Universal fault tolerant quantum computation on bilinear nearest neighbor arrays 
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Assuming an array that consists of two parallel lines of qubits and that permits only nearest 
neighbor interactions, we construct physical and logical circuitry to enable universal fault tolerant 
quantum computation under the [[7, 1, 3]] quantum code. A rigorous lower bound to the fault tolerant 
threshold for this array is determined in a number of physical settings. Adversarial memory errors, 
two-qubit gate errors, and readout errors are included in our analysis. In the setting where the 
physical memory failure rate is equal to one-tenth of the physical gate error rate, the physical 
readout error rate is equal to the physical gate error rate, and the duration of physical readout is 
ten times the duration of a physical gate, we obtain a lower bound to the asymptotic threshold of 
1.96 x 10" 6 . 

PACS numbers: 03.67.Lx 
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I. INTRODUCTION 

Quantum computation will potentially enable the effi- 
cient solution of computationally difficult problems such 
as factorization database searching [2j and quantum 
system simulation [3\. However, for any quantum com- 
puter system to be feasible it must permit the scalable 
implementation of fault tolerant quantum error correc- 
tion (FTEC) [IS, HQ- FTEC will exact a reduction in 
the rate of effective gate failure with each recursive level 
of encoding provided that the underlying physical failure 
rate is below some threshold [1, Q . The specific value of 
this threshold is influenced by the chosen quantum code, 
and the efficiency of circuits devised to implement this 
code, and by the spatial and operational limitations of 
the physical qubit array. 

Given the considerable interest in quantifying the 
experimental requirements associated with large-scale 
quantum computation, a great number of numeric and 
analytic threshold estimates have been undertaken [l(J. 
Typically, threshold estimates have permitted unlimited 
range interaction between arbitrary pairs of qubits. In re- 
ality, controlled long-range interactions of this type will 
be difficult to achieve. As many quantum computing pro- 
posals are based upon short-range interactions, a more 
valid assumption is that only nearest neighbor (NN) in- 
teractions are available [Tl|. The prevalence of linear 
nearest neighbor (LNN) proposals also suggests that a 
restriction may apply to the dimensionality of the qubit 
array in practice. 

Though a threshold has been shown to exist for a 
one-dimensional array with next nearest neighbor inter- 
actions [13], the viability of a LNN system remains an 
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open question. The primary difficulty associated with 
a strict LNN array is that a single fault in one of the 
SWAP gates used to interleave adjacent logical qubits 
prior to a transversal interaction can impart more than 
one error to a single logical qubit, violating the require- 
ments of FTEC under a distance-3 quantum code. One 
obvious remedy for this problem would be to implement 
a larger distance quantum code, such as the [[25, 1, 511 
Bacon-Shor code 1131; the [[23,1,7]] Golay code [IJ], or 
the [[11,1,5]] code jl5j |. which is the most compact known 
distance-5 code. In addition, the use of a LNN array 
does not preclude the implementation of a large quan- 
tum algorithm [lfij . However, due to the increase in the 
number of qubits and the complexity of the circuitry re- 
quired to detect and correct multiple-qubit errors, it is 
expected that such a restriction will incur a significant 
threshold penalty [ijj • Other ideas under investigation 
include a two-dimensional NN array, for which a favor- 
able fault tolerant threshold has been obtained 11811 • o r 
systems with a mechanism for qubit transport [19 . l2fj , 
though there is no clear path to the near-term fabrication 
of any such architecture. 

In this paper we consider a quasi one-dimensional array 
that consists of two parallel lines of qubits and permits 
only NN interactions. A bilinear architecture of this kind 
appears to be more tractable than any two-dimensional 
architecture and does not prohibit universal fault toler- 
ant quantum computation under a distance-3 quantum 
code. Such an architecture could be realized, for exam- 
ple,Jn a donor-based system, as was originally proposed 
in [2l|, or using superconducting technology [22]. We 
note that a CNOT threshold has been estimated for a NN 
array that increases in length and width with each level 
of concatenation [23| . In this paper we present physical 
and logical circuitry to achieve universal fault tolerant 
quantum computation on a more constrained bilinear ar- 
ray. In addition, we determine a rigorous lower bound to 
the threshold for universal fault tolerant computation on 
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this array in a number of physical settings. In the setting 
where the physical memory failure rate is negligible, we 
obtain a lower bound to the threshold that is an order of 
magnitude higher than that presented in [23]. 

The paper is organized as follows. In Section |TT] we 
review the chosen FTEC protocol and the method for 
threshold estimation. In Section HTTl we describe the adap- 
tation of non-local circuitry for the fault tolerant imple- 
mentation of a universal set of quantum gates to physi- 
cal and logical circuitry on a bilinear NN array. Section 
IIVI contains the calculation and presentation of thresh- 
old conditions for this array. Section fVl concludes with a 
summary of results and a description of further work. 



II. FAULT TOLERANT ERROR CORRECTION 

The seven qubit [[7, 1, 3]] code [24[ is sufficient to pro- 
tect against an arbitrary single qubit error and against 
multiple errors with non-zero probability of success and is 
an attractive code for concatenation as a number of fault 
tolerant logical gates can be implemented transversally. 
For example, as for any CSS code, a logical CNOT is re- 
alized by the parallel application of a CNOT from each 
qubit in the control logical block to each corresponding 
qubit in the target logical block. The set of logical sin- 
gle qubit gates {X, Z, H, S} can also be implemented in 
this fashion. This property leads to a reduction in the 
complexity of a number of fault tolerant gates relative 
to other CSS codes and to non-CSS codes where telepor- 
tation may be required. Furthermore, though resource 
requirements per logical qubit are slightly higher under 
the [[7,1,3]] code than under a block code [25|], under 
the [[7, 1, 3]] code it is possible to operate on each logical 
qubit simultaneously, potentially enabling faster compu- 
tation. 

The most gate- and time-efficient method of syndrome 
extraction under the [[7,1,3]] code is the encoded block 
method [26|. This method requires a seven qubit ancilla 
block to be prepared in the [[7, 1,3]] encoded state \0l)- 
The sequential application of transversal CNOT gates 
between the ancilla block and the data block allows the 
extraction of Z- and X-type syndrome information re- 
quired for FTEC. For example, the X-type syndrome is 
extracted by applying a transversal CNOT with the an- 
cilla block as the control and the data block as the target. 
This gate will copy any Z error on the data block to the 
corresponding location on the ancilla block. Measure- 
ment of the ancilla block followed by a classical parity 
check will determine the location of a single Z error. 

If the circuit for ancilla preparation is not fault toler- 
ant, a single fault during preparation of the state |0l) 
may lead to more than one error being copied to the 
data block. For example, during X-type syndrome ex- 
traction any X errors in the ancilla block will be copied 
to the data block. While a probabilistic procedure of 
ancilla verification is typically employed to restore fault 
tolerance, an alternative procedure has been proposed in 
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FIG. 1: X-type syndrome extraction using the encoded block 
method and ancilla decoding instead of verification. R is a 
recovery operation that may involve single qubit X and Z 
operations. 
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FIG. 2: To generate an encoded circuit each physical location 
(0-Ga) is replaced by a rectangle comprising the correspond- 
ing logical location (1-Ga) followed by error correction (1-EC) 
of each logical qubit. 
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FIG. 3: Two or more faults must occur within any extended 
rectangle for an encoded circuit to fail. 



which ancilla verification is replaced by a decoding circuit 
that is applied after the transversal CNOT between the 
data and ancilla blocks [27j ■ Though a single fault dur- 
ing ancilla preparation may still lead to a multiple-qubit 
error being copied to the data block, the decoding circuit 
is designed such that it is possible to determine the type 
and location of this error. A suitable recovery operation 
will take the form of a product of single qubit X and 
Z operations. FTEC using ancilla decoding is attractive 
as the qubits and circuitry that were required to verify 
the ancilla are no longer necessary. The circuit to im- 
plement X-type syndrome extraction using the encoded 
block method and ancilla decoding is shown in Fig. [TJ 

For an architecture that can interact arbitrary pairs 
of qubits, to generate an encoded quantum circuit each 
physical qubit is replaced by a logical qubit. Each lo- 
cation in the original circuit (gate, memory readout or 
state preparation) is also replaced by a rectangle [HI]. 
Each rectangle comprises a number of physical locations 
and executes a fault tolerant logical operation equivalent 
to the location type being replaced, followed by a FTEC 
cycle applied to each logical qubit involved in the loca- 
tion. This replacement procedure is shown schematically 
for a two qubit gate in Fig. [5J The resultant encoded 
circuit will have a higher effective fidelity than the orig- 
inal physical circuit provided that the failure rate of the 
physical locations is below some threshold. To increase 
the fidelity of the encoded circuit, this replacement pro- 
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FIG. 4: The single-error SWAP of computational qubits \cl) 
and |c2) on a vertical bilinear array. A physical SWAP gate 
is only ever applied between a computational qubit and a 
placeholder qubit. 
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FIG. 5: Logical circuit for ancilla encoding. Compound gates 
are used and a double-barred SWAP symbol indicates a single- 
error SWAP. 
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FIG. 6: Logical circuit for ancilla decoding. All measurements 
are made in the computational basis. 
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FIG. 7: Logical mesh circuit. The corresponding unmesh 
circuit is the circuit reversed. 



cedure can be applied recursively, whereby the new phys- 
ical locations are themselves replaced by rectangles. This 
recursive encoding forms a concatenated quantum code. 
Note that some technical details of this method need to 
be modified for the array considered in this paper. 

To estimate the fault tolerant threshold an extended 
rectangle is considered (28j . which is identical to a rect- 
angle but is preceded by a FTEC cycle applied to each 
logical qubit. An extended rectangle is shown schemati- 
cally in Fig. [3] Under the [[7, 1, 3]] code, assuming that all 
circuits are fault tolerant, a single fault in any extended 
rectangle will not lead to the failure of an encoded circuit 
as the resultant error or errors will be correctable. Two 
or more faults must occur within any extended rectangle 
for a circuit to fail. Assuming that every combination of 
more than one fault will lead to failure, a lower bound 
to the threshold is obtained by counting the number of 
each type of location within the circuit of each extended 
rectangle, as described in detail in Section [TV] 

III. PHYSICAL AND LOGICAL CIRCUITRY 

Adaptation of the chosen FTEC protocol to a bilinear 
NN qubit array involves the formulation of appropriate 
quantum circuits for ancilla encoding and for syndrome 
extraction. To achieve universal computation, circuits 



are also required for the fault tolerant implementation 
of a universal set of logical gates. One such set com- 
prises {I,H,X,Z,S,S\T,CNOT}. As the [[7,1,3]] code 
permits the transversal implementation of each of these 
gates except for T, and assuming that single and two 
qubit gates can be compounded, the circuits required for 
the error corrected universal set can be adequately rep- 
resented by {/, T, SWAP}, where / will be referred to 
as a memory location. This is possible because the cir- 
cuit for the SWAP extended rectangle comprises a similar 
(and slightly larger) gate layout as that for the CNOT 
extended rectangle. Therefore, to estimate the thresh- 
old for universal quantum computation, the construction 
of circuitry describing the extended rectangles of mem- 
ory, SWAP, T and readout is sufficient. Note that the 
readout extended rectangle includes transversal readout, 
qubit resetting and logical state preparation. 

In the absence of any other fault tolerant transport 
mechanism, for a NN architecture to be feasible it must 
include a fault tolerant SWAP mechanism at the phys- 
ical level and at each logical level. Under the [[7,1,3]] 
code, one fault during this process should impart no more 
than one error to any logical qubit. To satisfy this re- 
quirement, the qubit array is divided between computa- 
tional qubits, which include all data and ancilla qubits, 
and placeholder qubits, the states of which are inconse- 




FIG. 8: Logical circuit for the SWAP extended rectangle. Bold lines represent data qubits and regular lines represent ancilla qubits. Encode and decode components 
refer to Fig. [5] and Fig. [6] respectively. Triangular components are appropriately sized mesh and unmesh circuits similar to Fig. [7] Rectangular components represent 
various transversal interactions: 1 is a CNOT with each ancilla qubit the control and each data qubit the target; 2 is a CNOT with each data qubit the control and each 
ancilla qubit the target, and with H applied to each ancilla qubit before and after the CNOT; and 3 is a transversal SWAP between logical data blocks. Additional 
SWAP gates are included in transversal interactions where required. Vertical dots indicate the presence of qubits that are not required during the CNOT extended 
rectangle, but are required to achieve the fault tolerant T. The central mesh circuitry is not shown explicitly but all locations are included in the data presented in 
Table El 
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FIG. 9: Fault tolerant circuit for the T extended rectangle that is based on the implementation in Fig. 1101 EC is a full error 
correction cycle and ECx and ECz are partial error correction cycles that extract the X- and Z-type syndrome respectively, 
as shown in Fig.[T] Each CAT refers to the circuit in Fig. 1111 The preparation of the state \A*) = T\+) is as follows: ECx will 

prepare the state \0l) from the state |0), then two measurements of the operator SX — TXT^ are performed followed by a full 
error correction cycle. If an error is detected in the EC following CAT2 or in either of CATi or CAT2, or if the measurement 
results of CATi and CAT2 disagree, then the ancilla preparation is aborted, the ancilla is reset and the additional circuitry 
denoted by EXTi is used. To complete the preparation, logical Z is applied to the ancilla block conditional upon the eigenvalue 
of the operator SX = TXT^ being — 1. Following the interaction between the ancilla and data, logical S is applied to the 
data block conditional upon the associated measurement outcome being 1. If an error is detected by the second syndrome 
measurement within the first EC, and if EXTi has not already been applied, the additional circuitry denoted by EXT2 is used. 
If an error is detected in EXT2 the additional circuitry denoted by EXT3 is used. 
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FIG. 10: Implementation of the logical rotation T = 
exp( ~^ z ), where \Ax.) = T\+). Logical 5" is applied to the 
data block conditional upon the measurement outcome being 
1. The required state \Az.) is prepared by the circuit in Fig. 

eh 

quential to any computational processes 12]. Then, as 
shown in Fig. 2J the single-error SWAP of two compu- 
tational qubits is realized through a sequence of physical 
SWAP gates that move two computational qubits past 
each other on the bilinear array. During this sequence 
each physical SWAP is only ever applied between a com- 
putational qubit and a placeholder qubit. As the states 
of placeholder qubits are inconsequential, a single faulty 
SWAP will impart only one error to a single computa- 
tional qubit and, therefore, only one error to a single 
logical qubit. This forms what is effectively a SWAP- 
based fault tolerant transport mechanism. A single-error 
logical SWAP can be constructed from this underlying 
mechanism in order to satisfy the fault tolerant require- 
ment at all levels. At any level above the physical level, 
the bilinear NN physical qubit array becomes what is ef- 
fectively a LNN logical qubit array with the added prop- 
erty of single-error logical SWAP gates. Thus, in the 
formulation of circuitry describing the required extended 
rectangles, a distinction is made between logical circuitry 
and physical circuitry: 

Logical circuitry. With the introduction of single-error 
SWAP gates, LNN circuitry can be made fault toler- 
ant using a distance-3 code at any logical level. Log- 



FIG. 11: Circuit to prepare the required state \As.) = T|+). 
Logical Z is applied to the data block conditional upon the 
measured eigenvalue of the operator SX = TXT' being — 1, 
though in the circuit in Fig. [9] this is only after the eigen- 
value has been determined after two or three separate mea- 
surements. The cat state decoding procedure suggested in 
[27l ] is used to detect (but not correct) errors. Also note that 
a pair can be eliminated when this circuit is repeated in 
the circuit in Fig. [9] 

ical circuitry to implement FTEC simply follows from 
the equivalent non-local circuitry. Figures [5] and [6] show 
the [[7,1,3]] ancilla encoding and decoding circuits respec- 
tively where all circuitry is valid for any level above the 
physical level and where a single-error logical SWAP is 
indicated by a double-barred SWAP symbol. Additional 
LNN mesh and unmesh circuits, shown in Fig. [7l are 
used to interleave and then separate data and ancilla 
blocks prior to and after X- and Z-type syndrome ex- 
traction and to interleave adjacent logical blocks prior 
to any inter-logical transversal interaction. Circuitry for 
the SWAP extended rectangle is constructed from these 
basic components and is shown in Fig. [8] Construction 
of the memory, T and readout extended rectangles simi- 
larly follows from non-local circuitry. For the T extended 
rectangle, our starting point was the non-local fault tol- 
erant circuit presented in [28j which we then made more 
compact, reducing both the qubit and gate counts, the 
details of which can be found in Figs. [Will Finally, for 
the readout extended rectangle, the state \0l) is prepared 
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FIG. 12: An example of a physical mesh circuit that is more 
efficient than the circuit generated by simply replacing each 
SWAP in the LNN circuit in Fig. [7] with the circuit shown in 
Fig. [4] Four discrete steps are shown. 
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FIG. 13: Circuitry is arranged at all levels such that the du- 
ration of one T rectangle is equal to that of five memory 
rectangles, the duration of one readout rectangle is equal to 
that of two memory rectangles and the duration of one SWAP 
rectangle is equal to that of one memory rectangle. 



physical architecture. 



from the state |0) by the circuit in Fig. [TJ 

We note that a number of measurement outcomes will 
result in a larger T circuit, as shown in Fig.0 For exam- 
ple, an additional partial or full FTEC cycle is required 
prior to the classically controlled S. This is to prevent 
the combination of an X and a Z error on different qubits 
in a single logical block, which may arise due to a single 
fault during syndrome extraction, from transforming to a 
combination of a Y error and an X or Z error on different 
qubits. To be consistent with the rigorous assumptions 
made in Section IIV( and as our analysis requires that 
the depth of the T rectangle does not dependent on par- 
ticular measurement outcomes, we consider the largest 
circuit that can result from a single fault. 

Physical circuitry. To construct physical level circuitry 
from LNN circuitry it would be sufficient to replace ev- 
ery single-error logical SWAP by the single-error SWAP 
shown in Fig. 0] To obtain more gate- and time-efficient 
circuitry we note that the logical LNN mesh circuit has 
a much more efficient physical gate implementation, as 
shown in Fig. than that suggested by the naive im- 
plementation shown in Fig. 2] Additional efficiency was 
achieved by including gates between physical qubits that 
are oriented perpendicular to the direction of the bilin- 
ear array. This resulted in circuitry that is significantly 
more efficient than that presented in [23[. All physical 
level circuitry was optimized by hand in this manner. 

To facilitate analysis, at the physical level and at all 
logical levels the duration of rectangles describing all 
locations were arranged to permit parallel locations, as 
illustrated in Fig. 1131 Furthermore, since, in principle, a 
T could occur anywhere in a circuit at any time, during 
encoding every physical qubit is replaced by 21 physical 
qubits. It is expected that the detailed arrangement 
of gates and qubits will ultimately, however, be an 
automated process. This process may involve the dense 
packing of FTEC cycles during memory and may also 
incorporate routing to avoid any defective regions that 
are identified during some initial characterization of the 



IV. THRESHOLDS 

To obtain our threshold estimate it was assumed that 
faults would occur with a probability specific to the phys- 
ical location type and that this probability was equal for 
each location of the same type. It was also assumed that 
every combination of two or more faults would lead to 
a combination of errors that was not correctable under 
the [[7,1,3]] code. This assumption ensures a rigorous 
lower bound to the threshold. Note that it is possible to 
find and neglect pairs of faults that lead to correctable 
errors, thereby relaxing this pessimistic assumption and 
enabling a more precise determination of the threshold. 
This has been demonstrated in previous work to raise the 
lower bound by a factor of between 3.0 and 4.5 
Given these assumptions, the failure probability of any 
physical circuit is 

1 - (l-p m )^(l-p s )Ks {1 _ pr) N r 

- N mPm (l- Pm ) N <»- 1 (l-ps) IfB (l-Pr) If * m 

- (l-Pm^NsPsil-ps^-Hl-Pr)** lj 

- (1 -Pm) N ™(l -ps) Ns N rPr (l- Pr ) N r-\ 

where p m , ps and p r are the failure probabilities of physi- 
cal memory, gate and readout locations respectively, and 
where N m , N$ and N r are the numbers of physical mem- 
ory, gate and readout locations in the circuit. The sec- 
ond term corresponds to the failure of no locations. The 
third and fourth terms correspond to the failure of a sin- 
gle memory location and a single two-qubit gate location 
respectively, where applied single- and two-qubit errors 
are chosen adversarially. The final term corresponds to 
the failure of a single readout location, where a bit flip 
is applied to the classical result. Similarly, the failure 
probabilities of the physical level memory, SWAP, T and 
readout extended rectangles are given by 
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Plj = 1 - (l-P0m) Nomj (l-P0s) NoSj (l-P0r) Nor i 

- N 0mj p 0m (l ~ POrn)^- 1 ^ - PQS) N ° S] (1 - POr) N ° rj 

- (1 - POm) N °^N OSj p OS (l - posj^-^l ~ POr ) N ^ 

- {l- P a m ) N °^(l~POs) NoS] N 0r] por{l-pOr) N ^-\ 



where po m , Pas and Por are the probabilities of failure 
of physical memory, two-qubit gate and readout loca- 
tions respectively and where iVoy is the number of i-type 
locations in the physical level (level-1) j extended rect- 
angle for i = {m, S, r} (memory, SWAP, readout) and 
j = {m,S,T,r} (memory, SWAP, T, readout). To ac- 
count for the expected variation in timing of gates and 
readout, No m j is a expressed as a function of £ r , which is 
defined as the ratio of physical readout time to physical 
gate time. Therefore, as the relative duration of physical 



readout increases the number of memory locations in any 
extended rectangle will increase proportionately. 

The failure probabilities of logical level extended rect- 
angles can be expressed similarly, but as functions of the 
failure probabilities of the physical level extended rectan- 
gles, pi m , pis, Pit andpi r , as given by Eqs. (2-5), rather 
than of the physical locations. For example, the failure 
probability of the first logical level (level-2) T extended 
rectangle is 



PIT = 1 - (l-Plm) N ^ T (l-Pls) NlST (l-PlT)^ TT (l-pir) NlrT 

- N lmTPlm (l - Plm)^- 1 ^ ~ P1S) N1ST (1 - PlT) NlTT {l ~ Plr) N ^ T 

- (l-Pl m ) NlmT N 1S TPls(l-Pls) N ^- 1 0--PlT) IflTT (l-Plr) IfM (3) 

- (1 -pim) NlmT {l -P1S) N1ST N 1T TP1T(1 -Pit) Nitt ~ 1 {1 ~Plr) NlrT 

- {l~Pl m ) NlmT {l~Pl S ) NlST {l-piT) NlTT N lrT p lr {l~Plr) N ^ T -\ 



where Nut is the number of i-type locations in the logical 
level T extended rectangle for i = {m, S,T,r}. Polyno- 
mials that describe the logical level memory, SWAP and 
readout extended rectangles take a similar form to Eq. 
(6), though there are no T gates in any of these circuits. 
Due to the self-similarity of circuitry at all logical lev- 
els, polynomials describing level-n extended rectangles 
are identical to the corresponding level-2 polynomials, 
though they are expressed in terms of the failure proba- 
bilities of level- (n— 1) extended rectangles. Furthermore, 
as the of duration circuitry at all logical levels is artifi- 
cially set, there is no requirement for time-dependence 
to be included in any logical level polynomials. Tables U 
and [II] show values of iVoy and Nuj determined from the 
physical and logical circuitry presented in Section [IIII If 
memory locations are not neglected, the T extended rect- 
angle contains the most locations at the physical level and 
at all logical levels. As T is required to achieve universal- 
ity, the threshold for universal fault tolerant computation 
is given by the T threshold. 

Following the ideas presented in [29(, the level-n T 
threshold is where the probability of failure of the level-n 
T extended rectangle is equal to the probability of failure 
of the level-1 T extended rectangle. As any polynomial 
at any level can be expressed as a function of the fail- 
ure probabilities of physical locations by the recursive 
substitution of lower level polynomials, it is possible to 
determine a lower bound to the level-n threshold by solv- 
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TABLE I: Number of i-type locations in the physical level 
(level-1) j extended rectangle for i = {m,S,r} and j — 
{m, S, T, r}. The depth of each rectangle is also shown where 
the fundamental time scale is given by the duration of a phys- 
ical gate. 
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TABLE II: Number of i-type locations in the logical level-n j 
extended rectangle for i = {m, S,T,r} and j — {m, S,T,r}. 
The depth of each rectangle is also shown where the time scale 
is given by the duration of a level-(n— 1) memory rectangle. 

ing 

PnT (POS , POm , POr , U ) = PlT (POS , POrn ,P0r,U), (4) 

for pos, pom and po r . To estimate a physical gate thresh- 
old error rate it is necessary to specify the ratios of the 
physical memory failure rate and the physical readout 
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FIG. 14: p nT (pas, 0.1, 1.0, 10.0) for n = {1,2,3,4,5,100}. 
Lower bounds to the levels 2-5 and 100 T thresholds are 
1.36 x 10" 6 , 1.72 x 10" 6 , 1.85 x 10~ 6 , 1.91 x 10" 6 and 
1.96 x 10~ 6 respectively. 



failure rate to the failure rate of all other physical gates, 
denoted by R m and R r respectively. The specification 
of R m , R r and t r is referred to as a physical setting. 
Therefore, determination of a lower bound to the level- n 
threshold as an explicit gate failure rate in any physical 
setting involves solving 

PnT(P0S, Rm, Rr,t r ) = P1t(P0S , Rm, Rr , tr) , (5) 

for pos- Figure [T4l shows the failure rate of the levels 1-5 
and 100 T extended rectangles as a function of pos m 
the setting where R m = 0.1, Rr = 1-0 and t r = 10.0 and 
includes a lower bound to each of the levels 2-5 and 100 
T thresholds in this setting. Levels 2-3 thresholds are 
particularly useful as they indicate the maximum gate 
failure rates tolerable in a practical quantum computer. 
Note that if the gate failure rate is equal to the level-n 
threshold, if between 2 and n— 1 levels of encoding are 
used the encoded circuit will operate with a lower effec- 
tive fidelity than at only one level of encoding. Only 
with levels of encoding beyond n, and, therefore, only 
with greater resources, will a higher fidelity encoded cir- 
cuit be achieved. Realistically, gate failure rates must be 
at least an order of magnitude less than these thresholds 
to obtain significant benefit from FTEC. 

To make a comparison with thresholds previously ob- 
tained for other qubit arrays it is necessary to consider 
the asymptotic T threshold. For a physical gate failure 
rate equal to this threshold recursive encoding maintains 
the failure rate of the logical T extended rectangle. A 
lower bound to the asymptotic T threshold is approxi- 
mated by the level-100 T threshold, also included in Fig. 
[l~4l For a NN array that increases in length and width 
with each level of concatenation, in the setting where 
R m = 0.0, Rr = 1.0 and t r = 1.0 a CNOT threshold 
of 1.20 x 10~ 7 has been obtained [23[. In the identical 
setting, for a bilinear NN array we obtain a lower bound 
to the asymptotic T threshold of 2.88 x 10~ 6 . For a two- 
dimensional NN array, in the setting where R m = 0.1, 
R r = 1.0 and t r = 1.0, accounting for pairs of faults 



2.22x1 0" 6 
1.96x10 6 



POS 



0.10 

Rm. 



1.00 



FIG. 15: Asymptotic T threshold as a function of Rm in the 
setting where R r = 1.0 and t r = 10.0. 



Pos 




FIG. 16: Asymptotic T threshold as a function of R r in the 
setting where R m = 0.1 and t r = 10.0. 



Pos 




FIG. 17: Asymptotic T threshold as a function of t r in the 
setting where R m = 0.1 and R r = 1.0. 



that lead to correctable errors, a threshold for universal 
computation of 1.85 x 10~ 5 was obtained [l8j]. In the 
identical setting, for the bilinear NN array we obtain a 
lower bound to the asymptotic T threshold of 2.05 x 10~ 6 . 

Finally, to observe the dependency of the asymptotic 
threshold on the physical setting, R m , R r and t r were 
varied from the initial setting where R m — 0.1, R r = 
1.0 and t r = 10.0, as shown in Figs. [TMTTI respectively. 
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The greatest variation of the threshold is observed when 
R m is varied. This is expected as memory locations are 
the predominant type in most of the physical and logical 
level extended rectangles. However, the overall variation 
from the initial threshold is not great. For example, even 
assuming the extreme parameters R m = 1.0, R r = 100 
and t r — 1000, the threshold is reduced by less than two 
orders of magnitude to 3.78 x 10 -8 . 

V. CONCLUSION 

This paper presents physical and logical circuitry to 
achieve universal fault tolerant quantum computation on 
a bilinear NN array. The lower bound to the asymptotic 
T threshold derived from this circuitry is 1.96 x 10~ 6 . 
This result represents a significant improvement over 
a CNOT threshold previously obtained for a less con- 
strained array [231 ] - We note that this improvement is 
despite the requirement of universality and despite the 
application of a more rigorous method of threshold es- 
timation. As this improvement can be attributed to a 
reduction in the area of the physical and logical circuitry 
required to implement FTEC, it is expected that further 
advances in the efficiency of FTEC protocols and asso- 
ciated circuitry will correspond to further improvements 
in the threshold. 

While the results presented in this paper are highly rel- 



evant to existing NN architectures and to architectures 
that may be restricted to NN interactions between logi- 
cal qubits [U , further work is required to determine 
the threshold for universal fault tolerant computation on 
a LNN array under one or more larger distance codes. 
Such a result would help in assessing the viability of the 
large number of LNN architectures under development. 
Further work devising specific physical schemes of long- 
range interaction or transport that are fault tolerant is 
also desirable to permit the design of architectures with 
higher thresholds. 
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